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Abstract 

In this paper we show that the apparent failure of QCD lattice 
perturbation theory to account for Monte Carlo measurements of per- 
turbative quantities results from choosing the bare lattice coupling 
constant as the expansion parameter. Using instead "renormalized" 
coupling constants defined in terms of physical quantities, like the 
heavy-quark potential, greatly enhances the predictive power of lat- 
tice perturbation theory. The quality of these predictions is further 
enhanced by a method for automatically determining the coupling- 
constant scale most appropriate to a particular quantity. We present 
a mean-field analysis that explains the large renormalizations relating 
lattice quantities, like the coupling constant, to their continuum ana- 
logues. This suggests a new prescription for designing lattice operators 
that are more continuum-like than conventional operators. Finally, we 
provide evidence that the scaling of physical quantities is asymptotic 
or perturbative already at /3's as low as 5.7, provided the evolution 
from scale to scale is analyzed using renormalized perturbation theory. 
This result indicates that reliable simulations of (quenched) QCD are 
possible at these same low (3's. 
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1 Introduction 



In principle, nonperturbative lattice simulations allow the calculation of 
any quantity in QCD, without recourse to perturbation theory. In prac- 
tice, however, perturbation theory is important to lattice QCD in several 
ways. Firstly, it provides the essential connection between lattice simula- 
tions, which are most effective for low-energy phenomena, and the high- 
energy arena of perturbative QCD phenomenology. This is accomplished 
through such constructs as the operator-product expansion. Secondly, per- 
turbation theory can account for effects on low-energy phenomena due to the 
physics at distance scales shorter than the lattice spacing. Provided the lat- 
tice spacing a is small enough, systematic errors of order a and higher can be 
removed from the theory by perturbatively correcting the action and oper- 
ators that define the lattice theory. This approach provides a cost-effective 
alternative to simply reducing the lattice spacing when systematic errors 
must be removed. Finally, lattice simulations and perturbation theory must 
agree for short distance quantities, where both approaches should be reli- 
able, if we are to have confidence in simulation results for nonperturbative 
quantities. 

It is disturbing therefore that Monte Carlo estimates for most short- 
distance quantities seem to agree poorly with perturbative calculations. An 
example is the vacuum expectation value of the lattice gluon operator U in 
Landau gauge. This is the lattice analogue of the expectation value (A^) of 
the square of the bare gauge field A^. Since (Aj^) is quadratically divergent, 
the loop integral in first-order perturbation theory is dominated by momenta 
of order the cutoff, and perturbation theory should be effective for cutoffs 
of order a couple of GeV or larger. However, the perturbative result, when 
expressed in terms of the bare coupling constant oiat = di^t/'^'^ of the lattice 
theory, is 

(1 -iTr[/)pT = 0.97 aiat = . 078 (1) 



at (3 = = 6.llill This is almost a factor of two smaller than the 



nonperturbative result, 



(1 - \TtU)mc = 0.139, (2) 



obtained from Monte Carlo simulations.]^] The coupling constant is quite 
small here (aiat = 0.08), and the loop momenta typically large « vr/a 



(3 is the parameter used to specify the bare coupling constant in the standard lattice 
action for QCD. 
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6 GeV). Perturbation theory ought to work; instead it seems to fail com- 
pletely. Discouraging results such as this have arisen in a wide range of 
lattice calculations, leading to considerable pessimism about the viability of 
lattice perturbation theory at moderate /3's. 

In this paper we show that although these facts are true they are mis- 
leading. We find that the key problem with previous calculations of this 
sort is in the choice of the expansion parameter for the perturbation series: 
aiat is generally a very poor choice. There is no compelling reason in a field 
theory for using the bare coupling constant as the expansion parameter in 
weak-coupling perturbation theory. Standard practice is to express pertur- 
bation series in terms of some renormalized coupling constant, one usually 
defined in terms of a physical quantity. Indeed the renormalized coupling is 
usually a running coupling "constant" whose value in a particular expansion 
depends upon the length scales relevant in that process; there is no single 
expansion parameter for all series. The perturbative quantities important in 
lattice QCD generally involve lengths of order the lattice spacing a, and so 
one might expect little renormalization of the coupling from its bare value. 
However this argument, the usual rationale for using aiat, ignores the pos- 
sibility of a large scale-independent renormalization of the bare coupling. 
We find that just such a renormalization does occur in lattice QCD, making 
expansions in aiat useless except at very large /3's. 

Faced with large renormalizations, we must replace aiat by a renormal- 
ized coupling. It is straightforward to reexpress lattice perturbation ex- 
pansions in terms of any of the expansion parameters that have proven 
effective in continuum perturbation theory — for example, 0^3(9) with some 
physically motivated momentum q. When this is done, we find that lattice 
perturbation theory becomes far more reliable. In fact, perturbation the- 
ory becomes about as effective for lattice quantities as it is for continuum 



The large renormalization of aiat is due to the structure of the link 
operators from which the theory is built. The nonlinear relation between the 

'^'This situation in lattice theory parallels that for dimensional regularization, the other 
widely used regulator in QCD. Early calculations using dimensional regulatization were 
expressed in terms of the minimal-subtraction coupling constant ums, the "natural" def- 
inition for that regulator. The results usually looked nonsensical, with large coefficients 
appearing in the higher-order terms of most expansions. Consequently a modified minimal- 
subtraction scheme, the MS scheme, was introduced for defining the coupling constant. 
This scheme, while somewhat arbitrary, did result in reasonable perturbation series, and 
has since become standard. An analogous shift, away from aiat, is required in the study 
of lattice quantities. 
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link operator and the gauge field leads to large renormalizations of lattice 
operators relative to their continuum analogues, and these in term result 
in large shifts of the coupling constants in the action. In this paper we 
present a simple nonperturbative procedure for removing the bulk of these 
large "tadpole" renormalizations from gluon and quark operators. This 
procedure elucidates the problems with aiat • More importantly, perturbative 
expansions of the renormalization constants that relate quark currents and 
other composite operators on the lattice to their continuum counterparts 
become far more convergent once the tadpole contributions are removed. 

In Section 2 of this paper we discuss the symptoms that result from 
a poor choice of expansion parameter in a perturbation series. We show 
how these symptoms afflict lattice expansions expressed in terms of aiat, 
and we suggest a new, physically motivated procedure for renormalizing 
lattice perturbation theory.t^ll In Section 3, we compare predictions from 
our renormalized perturbation theory with nonperturbative results obtained 
from Monte Carlo simulations. We examine quark masses, {TrU), and a va- 
riety of Wilson loops and Creutz ratios. We find impressive agreement for 
all quantities, with no tuning of the theory, even at /3's as low as 5.7. In 
Section 4 we discuss the origins of the large renormalizations that arise when 
comparing lattice quantities with their continuum analogues. We develop a 
new prescription for building lattice operators that are much closer in behav- 
ior to their continuum counterparts; in particular the large renormalizations 



disappear.tlll The success of renormalized perturbation theory at low /3's 
suggests that the evolution of the coupling constant with lattice spacing is 
also perturbative and scaling asymptotic at these /5's. This important issue 
is discussed in Section 5. Finally, in Section 6, we summarize our conclu- 
sions, stressing their implications concerning the reliability of simulations 
on relatively coarse (and therefore much less costly) lattices. The data for 
the plots throughout the paper is tabulated in the Appendix. 

2 Renormalized Lattice Perturbation Theory 

2.1 A poor expansion parameter 

If an expansion parameter Ogood produces well behaved perturbation se- 
ries for a variety of quantities, using an alternative expansion parameter 



(3)pQj. discussion of these issues in the context of dimensionally regularized QCD per- 
turbation theory, see Q|. A preliminary version of our lattice results is in 
'^'A preliminary version of this analysis is published in El. 
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Obad = agood(l " 10,000agood) will lead to second-order corrections that 
are uniformly large, each roughly equal to 10,000abad times the first order 
contribution. Series expressed in terms of abad; although formally correct, 
are misleading if truncated and compared with data. The signal for a poor 
choice of expansion parameter is the presence in a variety of calculations 
of large second-order coefficients that are all roughly equal relative to first 
order. 

A large coefficient appears in the first second-order calculation done on 
the lattice: the calculation of the gluonic three-point function used to relate 
the A parameter of the bare lattice coupling aiat to the A's of various con- 
tinuum coupling constants.!^, ^ The coupling constant Ci{q)^;^^, defined in 
terms of this three point function at momentum q, has the expansion 

«(9)r55Si = «iat {l + aiat(/3o H^^/aqf + 5.419)} , (3) 

where /3o = ll/47r. Naively, one expects that a{q = Tr/a)^^^ « aiat, since 
vr/a is roughly the largest momentum on the lattice. The constant 5.419 
spoils the equality; it results in very large ratios between continuum and 
lattice A's. 

Since continuum quantities are usually well behaved when expanded in 
terms of it is immediately obvious that most other continuum 

quantities will have a similar constant term when expressed in terms of aiat • 
For example, the heavy quark potential V{q) at momentum transfer q has 
the expansion 0] 

V(,) = {l + a,., (a In (i)' + 4.70) } (4) 

where Cf = 4/3 is the quark's color (Casimir) charge. Similar results hold 
for the e^e~ hadronic cross section, derivatives of moments for deep inelastic 
ep scattering, etc. 

A crucial point is that a similar constant term appears in the expansions 
for all short-distance lattice quantities that have been studied. For example, 
the corrections to the heavy-quark potential as a function of distance have 
the form|p 

V{R) = |l + aiat (^/3oln (^)' + C{R/a)^ | , (5) 

where C{R/a) for various values of R is given in the following table: 
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R/a 


2 4 


6 oo 


C{R/a) 


5.5 5.5 


5.6 5.711 



(The constant for R = oo can be obtained by Fourier transforming the 
equation for V{q) above.) Note that the constants C{R/a) at finite R vary 
httle from the one at R = oo. This is expected since these corrections are 
dominated by quadraticahy UV divergent tadpole loops that are insensitive 
to the external momenta. 

As we show later, similar terms are present in Wilson loops and Creutz 
ratios. Thus the pattern of second-order coefficients for lattice quantities 
strongly suggests that aiat is a poor choice of expansion parameter. 



2.2 A better expansion parameter 

To define an improved (renormalized) expansion parameter, we must both 
choose a definition of the running coupling as{q) ("fix the scheme") and 
specify how the scale q of the coupling is to be chosen ("set the scale"). 
It is natural and convenient in perturbation theory to tie the scale of the 
coupling to that of the loop momenta circulating in the Feynman diagrams. 
Thus, we want to define as{q)_sp that it approximates the coupling strength 
of a gluon with momentum q!(^ It is also important that as{q) be defined in 
terms of a physical quantity, so as to avoid confusions, such as that between 
the MS and MS schemes, that are artifacts of arbitrary definitions. 



2.2.1 Fixing the Scheme 

Of the many physical quantities one might use to define an as{q), the heavy- 
quark potential V{q) is among the most attractive. |§] Typically there is an 
integral over the momentum of the leading-order gluon, but the gluon in 
V{q) has only momentum q. Thus it is particularly easy to tie the coupling 
constant's argument to the gluon's momentum for this quantity: we define 
av{q), the coupling strength of a gluon with momentum q, such that 

Viq) = (6) 



^^^It is natural in a gauge theory to associate the scale of the coupling with the gluon's 
momentum since every g in the theory is associated with a particular by gauge invari- 
ance. This association allows us to set the scale in a gauge invariant way. 
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with no higher-order corrections. We can easily relate ay to the bare lattice 
coupling constant aiat since V^q) has been computed in terms of aiat (Eq. 

aiat = OiY{q) {l - ay (/5o ln(7r/ag)2 + 4.702) } + 0{(x\). (7) 

for SU(3) color with no light-quark vacuum polarization. With this expres- 
sion, any one-loop or two-loop lattice expansion can be reexpressed as a 
series in ay. The q dependence of ay is given by the usual formula, 

ay\(l) = Po Hq/^v? + /3i//3o l^Hq/Kyf + 0{ay{q)), (8) 

where I3q = ll/in (as before), /3i = 102/167r^, and 

Ay = 46.08 Aiat (9) 

is the scale parameter for this scheme. The scale parameter Aiat for oiat is 
defined implicitly by 

a^^l = Po ln(l/aAiat)' + Pi/Po lnln(l/aAiat)' + • • • ■ (10) 

Note that ct^Asiq) is numerically fairly close to ay{q), and thus is another 
useful alternative to aiat- In this case, 

aiat = aMs(g) {l - «MS (/5o ln(vr/ag)2 + 3.88o) } + 0(a|g), (11) 
and the scale parameter is 

Aj;jg = 28.81 Aiat . (12) 

2.2.2 Setting the Scale 

The coupling constant ay is defined so that ay{q*) is the appropriate ex- 
pansion parameter for a process in which the typical gluon momentum is q*. 
For many processes it is possible to guess q* fairly accurately. For example, 
power-law UV divergent quantities like {TiU) are controlled by the lattice 
modes with the highest momenta, and so one expects q* ~ n/a. Although 
such guesses are often sufficient, there is a simple automatic procedure that 
takes the guessing out of q*. Such a procedure has proven invaluable in our 
systematic study of the reliability of perturbation theory. 

Consider a one-loop perturbative contribution in our scheme: 

I = ayiq*) J d\f{q) (13) 
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where q is the gluon's momentum. The natural definition of q* would be 



av{q*) j d\f{q) = j d\av{q)f{q) (14) 

except that the second integral is singular. The singularity is due to the 
pole in the coupling constant at g = Ay. This pole is an artifact of the 
all orders summation of perturbative logarithms that is implicit in the for- 
mula for av{q) (Eq. ^). The singularity does not arise in any finite order 
of perturbation theory, as may be seen by replacing the running coupling 



constant av(^) in Eq. 14 by its expansion in terms of the coupling constant 



renormalized at some fixed scale /i: 

av{q) = o^v{^l){l + Po ln{q / fif ay (fJ^) + (/?o ln(g///)^av(//))^ H } (15) 

None of these terms separately results in a singularity, but the sum of all 
terms diverges. 

In fact it is incorrect to sum to all orders since the QCD perturbation 
series is an asymptotic series. The proper procedure is to retain only those 
terms consistent with the accuracy of the rest of the calculation. For our 



purposes we should retain only the first two terms in Eq. 15 
av{q*) I d\f{q) ^ ayif^) f d\f{q) + 



Poavif^rj d\f{q) \n{q/f,Y + ---. (16) 

Expanding av{q*) in terms of ayifJ-) in this equation, we obtain a simple 
definition for q* (independent of fi): 

ln(«*2) = /^VMM^. (17) 

2.2.3 Summary 

To summarize, our general procedure for analyzing a perturbation series in 
lattice QCD involves replacing aiat by av{q*) using Eq. |^ The scale q* is 
determined by probing the first-order calculation with a factor Ing^, as in 



Eq. [17|. In calculations that extend through two-loop order, we assume that 
the one- loop q*, determined this way, is also appropriate for the two-loop 
contribution. 

A special feature of expansions in av{q*) is that they are unaffected 
through second order by quark vacuum-polarization insertions in the gluon 
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propagator. All such contributions are automatically absorbed into av{q*), 
by virtue of its definition. As a consequence all of the perturbative ex- 
pansions we use in this paper (second order as well as first) are identical in 
quenched and unquenched versions of QCD when they are expressed in terms 
of ay(g*). Only the evolution of av{q) is different: Pq — > (11 — |nj)/47r 
and Pi — > (102 — 18nj)/167r^ in Eq. |8| for nf light-quark flavors. The Uf- 
independence of oy expansions leads to an alternative procedure for deter- 
mining q* that is analyzed extensively, for continuum QCD, in ||2|. 

3 Testing Renormalized Perturbation Theory 

Our procedure for defining a renormalized coupling constant with a proper 
scale (Section p.2| ) follows solely from known results in lattice perturbation 
theory, without regard to Monte Carlo data. Only now are we ready to 
consider the extent to which our renormalized perturbation theory agrees 
with Monte Carlo simulations of short-distance quantities. 

Having converted all of our perturbative expansions from ctiat to ay, we 
need some way of determining values of ay that are appropriate to particular 
simulations. The most straightforward procedure is to measure ay in the 
simulationsE^ This can be done, for example, by measuring the heavy-quark 
potential, or, more simply, by measuring the trace of the plaquette operator 
Upiaq (the 1x1 Wilson loop). The improved perturbative expansion for the 
logarithm of Trf/piaq is 

- ln(iTrC/piaq) = 4.18879ay(3.41/a) {l - l.lDay + ©(ay)} . (18) 

Given data for this quantity, one can easily solve for ay(3.14/a). The cou- 
pling ay{q) for other g's can then be obtained using standard two-loop 
evolution (Eq. |8|). We have extracted o;y(3.41/a) in this way from data for 
quenched QCD at several /3's. The results, evolved down to q = 1/a, are 
given m Table |l[ We also give values for aiat and for aj^{l/a), the latter 
being obtained from the measured ay using the relation = 0.6252 Ay. 

Our choice of — ln(iTrC/piaq) for determining ay is for convenience; 
we have not attempted to optimize this choice. One could use any other 
short-distance quantity whose ay-expansion is known through second order. 

'^'As we discuss in Section 4, ay can also be computed directly from aiat without using 
Monte Carlo data. However this procedure is probably less accurate than measuring ay, 
particularly at lower /3's. 
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/3 


-InPFii 


aiat 




av(l/a) 


5.7 


0.5995 


0.0838 


0.2579 


0.3552 


6 


0.5214 


0.0796 


0.1981 


0.2467 


6.1 


0.5025 


0.0783 


0.1860 


0.2275 


6.2 


0.4884 


0.0770 


0.1774 


0.2144 


6.3 


0.4740 


0.0758 


0.1690 


0.2020 


6.4 


0.4610 


0.0746 


0.1617 


0.1913 


9 


0.2795 


0.0531 


0.0815 


0.0878 


12 


0.1954 


0.0398 


0.0532 


0.0558 


18 


0.1227 


0.0265 


0.0317 


0.0326 



Table 1: Monte Carlo data for logarithm of the plaquette, together with the 
couphng constant values used in this study 

Other alternatives might be Creutz ratios of small loops, whose perturba- 
tive expansions might be more convergent, or a combination of Wilson loops 
chosen so that potential nonperturbative area-law contributions cancel (eg, 
lnTr[/2x2 -41nTr;7ixi). 

Note that, as discussed earlier, the formula used in measuring ay (Eq. |l^) 
is valid also for unquenched QCD, as are all of the ay expansions that follow. 
Thus precisely the same techniques and tests we use here can be applied to 
the unquenched case. We have not yet done this, but we expect similar 
results. 

3.1 {AD 

The lattice equivalent of (Aj^) is (1 — ^TrU), which is given in perturbation 
theory by 0.97 a^. The one-loop contribution comes from a quadratically 
divergent tadpole graph, and we therefore expect that it is dominated by 
momenta of order the lattice cutoff vr/a. Using the procedure of Section 
we find q* = 2.80/a. In Fig. |l| we compare perturbative results results for 
(1 — ^TrC/) with Monte Carlo dataQ at several values of /9. We present 
results from perturbation expansions in aiatj in our favorite coupling con- 
stant av{q*), and in ay^{q*). The data agree with perturbation theory to 
within 10-15% for all f3 > 5.7 when ay or is used. Uncalculated terms 
of order ay or higher in the perturbation theory could easily account for the 
remaining differences; the differences between the ay and a^jg predictions 
give an indication of how important such terms might be. Of course part 
of the difference between perturbation theory and Monte Carlo might be 
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Figure 1: The expectation value of the trace of a hnk in Landau gauge, 
calculated by Monte Carlo (circles), and in first order perturbation theory 
using for the expansion parameter av{q*) (diamonds), ay^{q*) (boxes), and 
aiat (crosses). 

nonperturbative, particularly at the lowest /?. The data disagree with the 
aiat exapnsion by almost a factor of two. 



3.2 Mass renormalization for Wilson quarks 

A famous example of the "failure" of lattice perturbation theory is the cal- 
culation of the additive mass renormalization for Wilson quarks. The bare 
mass in Wilson's formulation of the lattice quark action is renormalized by 
an additive power-law divergent term. The critical quark mass, for which 
this term is canceled (leaving the quark massless), is given in perturbation 
theory by rrica = 1/2kc — 4 = — 5.457as.||9| (Here, k is the "hopping param- 
eter" used to parameterize the quark mass in lattice gauge theory.) The 
linear divergence in this result suggests that the important momenta here 
are of order vr/a. We find q* = 2.58/a using our procedure (Eq. |^). In 
Fig. § we compare perturbative results for rric with Monte Carlo data|lC] 
at several values of /3. Again we see that the data agree with our renor- 
malized perturbation theory to within 10-15% for all /3's, but disagree with 
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Figure 2: The critical quark mass rric for Wilson quarks, calculated by Monte 
Carlo (circles), and in first order perturbation theory using for the expansion 
parameter av{q*) (diamonds), Oi-^{q*) (boxes), and aiat (crosses). Statis- 
tical errors in the Monte Carlo results are negligible. 

perturbation theory using Qiat by almost a factor of two. 



3.3 Wilson loops and Creutz ratios 

Aside from the heavy-quark potential and the coupling constant, Wilson 
loops are the only lattice quantities for which two-loop perturbation theory 
has been calculated. Consequently they provide the most stringent tests of 
perturbation theory. Large Wilson loops have badly behaved perturbative 
expansions for a trivial reason: they contain a self-energy contribution pro- 
portional to the length of the loop. For large loops, contributions to this 
self-energy approximately exponentiate, so we expect that the logarithm o f a 
Wilson loop is better behaved in perturbation theory than the loop itself. 



(7) 



^^'Our data confirms that perturbation theory works better for logarithms of the Wmn 
than for the Wmn themselves, the expansions for the latter failing completely for even 
modestly large loops. Curiously the pathologies in the Wmn expansions seem to cancel 
the pathologies in aiat when m and n are small, making the oiat expansion more accurate 
than the av expansion for these loops. Neither expansion is as accurate as expanding 
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Taking Creutz ratios[ll| of Wilson loops further improves perturbation the- 
ory by reducing the effects of both the divergent contributions associated 
with the perimeter of the loop and those coming from the corners of the 
loop. 

For these reasons, we concentrate in this study on the logarithms of small 
Wilson loops and on Creutz ratios Xmn defined by 

X^n = -\nl— — . (19) 

where Wmn is one third the expectation value of the trace of the m x n 
planar Wilson loop: 

= liTlUmxn). (20) 

We compare perturbative predictions with new data generated on a 16^ 
lattice at /3's ranging from 5.7 to 18.[|l2| We use one-loop and two-loop 
perturbation-theory coefficients computed for a 16^ lattice^], and include 



the leading-order contribution from the zero mode. [ 13 1 Thus our perturba- 



tion theory is accurate up to uncalculated terms of order ay, and of order 
Oy/V, due to the zero mode, where V = 16^ is the volume of the lattice. 
The finite-volume errors becomes significant for larger loops and so we limit 
ourselves to 5 x 5 loops and smaller. 

In Fig. ^ we show results for X22 , calculated through first order in ag , and 
also through second order. The pattern at first order is similar to that in our 
previous examples: expansions in av{q*) and a\a,tiQ*) give reliable results at 
all /3's; the expansion in aiat is off by almost a factor of four at /? = 5.7, and 
still by almost 30% at /? = 12. The second-order corrections significantly 
improve agreement between the data and the ay and expansions, with 
errors ranging from a few percent at /3 = 5.7 to a few tenths of a percent 
at f3 = 12. The remaining discrepancy could easily be accounted for by 
uncalculated corrections of order a^, although again nonperturbative effects 
may play a role at the lowest /?'s. The second-order expansion with aiat 
gives results that are at least an order of magnitude worse than those from 
the other two expansions (at all /3's). By comparison with the others, the 
convergence of this expansion is very sluggish — a unambiguous symptom of 
a bad expansion parameter.!!^ 



— In Wmn in powers of av and then exponentiating. This last procedure gives good results 
(when f3 is large) for all loops out to 8 x 8, the largest we checked. 

'*'Note that some of our results have been anticipated in the literature. The fact that 
perturbative results for Creutz ratios are better behaved when expanded in terms of ctj^ 
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Figure 3: Results for Crcutz ratio X22 at different couplings [3 from Monte 
Carlo simulations (circles), and from perturbation theory (using ayiq*) (di- 
amonds), oi^{q*) (boxes), and aiat (crosses)). The first plot shows per- 
turbation theory through one- loop order, and the second through two-loop 
order. Statistical errors in the Monte Carlo results are negligible. 
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Figure 4: Results from perturbation theory (with <y.Y(q*) (diamonds), 
aj^(g*) (boxes), and aiat (crosses)) and Monte Carlo simulations (circles) 
for diagonal Creutz ratios Xn,n at /? = 6.2. Statistical errors in the Monte 
Carlo results are negligible. 
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Figure 5: Results for the Creutz ratio X22 from Monte Carlo simulations 
(line) and from the perturbation expansion in av{q) (diamonds) versus the 
scale q. 

In Fig. ^ we show two-loop results with each of the coupling constants for 
a variety of different Creutz ratios at /5 = 6.2. The ay and expansions 
are again far superior for all of the ratios. 

We expect smaller momentum scales for Creutz ratios than for the loops 
themselves since many of the divergent contributions to loop expectation 
values cancel in the ratios. Our scale setting procedure indicates that q* is 
1.09/a for X22, and smaller for ratios involving larger loops. 

The importance of choosing a proper q* is illustrated in Fig. ||, where 
the two- loop prediction for X22 has been reexpressed in terms of ay(g) and 
plotted versus q. Taking q = vr/a, for example, rather than q = q* = 1.09/a 
results in a 10% error rather than a 1% error. This situation should be con- 
trasted with that for — lnTy22 (Fig. P). This quantity is significantly more 
ultraviolet than X22) having q* = 2.65/a. Here our perturbative estimate 
degrades significantly if we use, say, q = 1/a rather than q* to set the scale of 

than when expanded in terms of aiat was pointed out in . The fact that perturbative 
results for Creutz ratios are better behaved when expanded in terms of an defined from 
any given ratio than when expanded in terms of aiat was pointed out in jlSj. 
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In W22 




Figure 6: Results for — In W22 from Monte Carlo simulations (line) and from 
the perturbation expansion in av{q) (diamonds) versus the scale q. 

our expansion parameter. These examples illustrate the importance of our 
scale-setting procedure when high precision is required. Small departures 
from q* are unimportant, at least for reasonably convergent series; but de- 
viations by factors of two or more can affect the reliability of a perturbative 
estimate. 

4 Mean field theory 

We have shown that perturbation theory works well when a proper coupling 

constant is used, but it is still important to understand the origins of the 
large mismatch between the lattice coupling and the continuum couplings. 
This mismatch is one of many examples where a large renormalization is 
required to relate a lattice quantity to its continuum analogue. In this 
section we explore the connection between operators on the lattice and in 
the continuum. 
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4.1 Tadpole improvement 



We usually design lattice operators by mapping them onto analogous oper- 
ators in the continuum theory. For gauge fields, this mapping is based upon 
the expansion 

U^{x) ^ e^'^s^'^(^) ^ 1 + iagA^ix). (21) 

This expansion seems plausible when the lattice spacing a is small, but 
it is misleading since further corrections do not vanish as powers of a in 
the quantum theory. Higher-order terms in the expansion of f/^ contain 
additional factors of gaA^, and the A^'s, if contracted with each other, 
generate ultraviolet divergences that precisely cancel the additional powers 
of a. Consequently these terms are suppressed only by powers of g'^ (not 
a), and turn out to be uncomfortably large. These are the QCD tadpole 
contributions. 

The tadpoles spoil our intuition about the connection between lattice op- 
erators and the continuum, and so we should not be surprised if the lattice 
theory is not quite what we expected (because of large renormalizations) . 
In order to regain this intuition we must refine the naive formula that con- 
nects the lattice operator to the continuum operator (Eq. |2l| ). Consider the 
vacuum expectation values of these operators. In the continuum, the expec- 
tation value of 1 + iagAfj_{x) is 1. In the lattice theory, tadpole corrections 
renormalize the link operator so that its vacuum expectation value (in, say. 
Landau gauge) is considerably smaller than 1 (see Fig. |l|). This suggests 
that the appropriate connection with continuum fields is more like 

U^{x) ^uo{l + iagA^{x)), (22) 

where uq, a number less than one, represents the mean value of the link. 
Gauge invariance requires that parameter uq enter as an overall constant .li^ll 

'■^^This formula follows simply from a renormalization-group argument. The tadpole 
contributions come mainly from the gauge-field modes with the highest momenta. Con- 
sequently the tadpoles can be removed by splitting the gauge field into ultraviolet (UV) 
and infrared (IR) parts (in a smooth gauge), and integrating out the UV parts. Averaging 
over the UV modes, the link operator is replaced by its IR part: 

ito e*"'''^^' ^ Ui,{l + iagA]^), 

where now the Taylor expansion of the exponential is quite convergent. Parameter uq 
contains the averaged UV contribution. It enters only as an overall constant since the link 
operator functions as a gauge connection both before and after averaging. 
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The mean-field parameter uq depends upon the parameters of the the- 
ory. It can be measured easily in a simulation. Simply measuring the link 
expectation value gives zero since the link operators are gauge dependent. 



However relations such as Eq. 22 only make sense in smooth gauges, like 
Landau gauge. Thus one might define uq to be the expectation value of the 
link operator in Landau gauge. A simpler, gauge-invariant definition uses 
the measured value of the plaquette: 

uo = (iTrf/piaq)^/^ (23) 

Several other definitions are possible; all give similar results. At /? = 6, for 
example, uq is 0.86 from the Landau-gauge link and 0.88 from the plaquette. 



Our improved relation, Eq. |22|, between lattice and continuum gauge-field 
operators suggests that all links that appear in lattice operators should 
be replaced by C/^/tto, where uq is measured in the simulation. The opera- 
tors U^/uq are much closer in their behavior to their continuum analogues; 
large tadpole renormalizations are largely canceled out by the uq (and the 
cancellation is nonperturbative since uq is measured rather than calculated) . 
This is the key ingredient in our tadpole-improvement procedure for lattice 
operators. Several illustrations follow in succeeding sections. 



4.2 ay from aiat 

Our new prescription for building continuum-like operators suggests that 

4luon = E :F^'I^(f^plaq + h-C). (24) 

is a better gluon action for lattice QCD. In particular, perturbation theory 
in is much more like continuum perturbation theory (ie, no tadpoles). Of 
course this tadpole-improved action becomes the normal gluon action if we 
identify 

= glJ{\T,{U^,^)). (25) 

This is a very important relationship; it tells us that the correct expansion 
parameter for the usual theory is g^ rather than gf^^^. The difference is 
significant: for example, g^ ~ l.lgf^^ at /? = 6 (using the measured value 
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of the plaquette to relate the mnplings V T t is a big mistake to expand in 



(10) 



powers of aiat rather than diat = 5^/47r. 

If our mean- field analysis is correct, diat should be roughly equal to 
ai/(7r/a). This is confirmed by perturbation theory which implies that 

ay(vr/«) = n^u \ + 0-513ay + (26) 

(§TrC/piaq) 

the difference between the two coupling constants is only a few percent at 
P = 6. This formula provides a nonperturbative relationship between the 
bare lattice coupling aiat and ay when measured values for the plaquette 
are used. 

Note that since the renormalization is multiplicative, its main effect is to 
rescale the argument of the running coupling constant. This suggests that 
we define ay by 

av(46.08/a) = f 1 + 0{al)) (27) 



(since Ay = 46.08Aiat), and then use two-loop evolution to determine avig) 
for any other scale q. This provides a purely perturbative relation between 
ay and aiat- 

In Fig. 1^ we compare the measured values of ay(7r/a) from Section ^ 
with values obtained from the mean-field formula, Eq. p^, and from the 



perturbative formula, Eq. Large coupling-constant renormalizations are 
automatically incorporated when ay is measured, and so the validity of our 
mean-field analysis is tested by the extent to which the mean-field values 
agree with the measured values. All three methods produce results consis- 
tent up to corrections of order ay. The first two methods of determining 
ay are probably preferable to the perturbative formula at low /3's since they 
incorporate some higher-order and nonperturbative effects. 

Our prescription for defining tadpole- improved lattice operators is crucial 
in other, related contexts. One example is in defining operators to represent 
the chromoelectric and magnetic fields. These are needed for the operators 



'^°'Our analysis of the gauge-field action was anticipated by Parisi who gives a simple 
analysis for the compact abelian theory. To see what effect the (UV- divergent) tadpoles 
have on infrared modes, we can split the gauge field in UV and IR components, and 
average over the UV part. Then the abelian gauge action becomes 

{g-^ cos(5F,r + 5i^™))"'' = g'^cosigF^,'')) cos(3F™) 

and the effective coupling for the IR modes is g^ divided by the UV part of the plaquette 
expectation value. 
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Figure 7: Values of Q;y(7r/a) as determined by measuring — ln(iTr?7piaq) 
(circles), by using a nonperturbative mean- field formula to relate it to the 
bare coupling (diamonds), and by using perturbation theory to relate it to 
the bare coupling (crosses). 
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that remove 0{a,a^) errors from quark actions. The standard cloverleaf 
definitions, Ed and Bd, involve a products of four link operators, just like 
the plaquette [^. Thus the tadpole-improved operators, 

Eel = Eel/(iTV(C/piaq)) (28) 
Bel = Bel/(iTr(C/piaq)), (29) 

are almost twice as large at /3 = 6. The plaquette factors account for the bulk 
of the very large renormalizations found in perturbation theory for operators 
containing cloverleaf fields. Such operators play an important role in all 
formulations of heavy-quark dynamics; omitting the tadpole renormalization 
leads to severe underestimates of their effects. 



4.3 Improved Wilson fermions 

Our tadpole-improvement scheme provides valuable insights into the pat- 
tern of large renormalizations in lattice QCD, and it is generally trivial to 
implement. As another example consider the tadpole- improved action for 
Wilson quarks: 

5, = ^V^^ + K^V^('(l + 7^)^)^ + .... (30) 

Again, this action is identical to the usual one if we relate the modified 
parameters, here the hopping parameter k, to the usual ones by rescaling 
with uq: 

R = kuq. (31) 

The modified hopping parameter should be more continuum-like; for 
example, the tree- level value that gives massless quarks, kc = 1/8, should 
be roughly correct for interacting quarks as well, at least at high /3's. Thus 
an approximate nonperturbative formula for the critical value of the usual 
hopping parameter is 

Kc^l/SuQ. (32) 

This formula accounts for about 75% of the renormalization of the hopping 
parameter when /? is large, as is evident if we compare the perturbative ex- 
pansions for the two sides. By combining these perturbative expansions, we 
obtain a tadpole-improved perturbation theory for the critical bare mass nic 
(where rrica = 1/2kc — 4, as in Section 3.2): 

= -4 (l - (iTrC/piaq)^/^) - 1.268 ay (1.03/a) + ©(a^). (33) 
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/9 


Srric a 


perturbative 


tadpole- improved 


measured 


5.7 


-0.44 


-1.11 


-1.00 


-1.04(2) 


6.0 


-0.31 


-0.91 


-0.80 


-0.80(2) 


6.1 


-0.29 


-0.86 


-0.76 


-0.78(2) 


6.3 


-0.22 


-0.79 


-0.70 


-0.70(2) 



Table 2: Mass renormalization for Wilson fermions at different couplings /3 
as computed using ordinary ay perturbation theory, tadpole-improved 
perturbation theory, and Monte Carlo simulation. Also listed is the 
first-order perturbative correction to the mean- field estimate of rric, 
5mca = -1.268Qv(1.03/a). 



Using the measured value of the plaquette operator, this formula should 
be more accurate than the purely perturbative formula used in Section 
since large tadpole renormalizations are being summed to all orders. Higher- 
order perturbative corrections should be smaller for the improved formula, 
as should nonperturbative effects. This seems to be the case as we show in 
Table ^, where the two predictions are compared with Monte Carlo data. 
Tadpole-improved one-loop perturbation theory predicts rric (and Kc) about 
as accurately as it can be measured. 

In the continuum limit, the tadpole-improved lagrangian for massless 
quarks becomes 

2klp-ff,d^'i^ + 0{a). (34) 

This indicates that is the lattice quantity that corresponds to the 

continuum quark field. Since k ~ 1/8 when the quarks are ma ssless. a 
tadpole-improved operator for massless quarks on the lattice ii^^^t 

= ip/l. (35) 

This lattice operator has roughly the same normalization as the continuum 
field; in particular, there are no large tadpole contributions to the renormal- 
ization constant relating them. This is important in designing new lattice 
operators involving quark fields. For example, if one wants matrix elements 



(ii)'pjjjg result is only valid for k « Kc- Away from Kc, ip becomes Q 

l/i EE (1 - 3/t/4Kc)^''^ tp. 

Note that the k dependence here is quite different from that of the commonly used (but 
incorrect) 
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of the continuum current ■07^7^ -0) then one should simulate with the lattice 
operator 

^^M^5^ ^ l^^M^5^_ (36) 

Our procedure differs somewhat from common practice. Frequently the 
quark renormahzation factor is taken to be for massless quarks, rather 
than the factor \/2kc = 1/2 that we use. The former factor differs signifi- 
cantly from 1/2 unless /? is quite large. In our mean-field analysis, the con- 
ventional factor y/2nc is in effect divided by ^/8k^. This additional factor 
removes the bulk of the large tadpole corrections usually found in calcula- 
tions of renormahzation constants for quark operators. We have verified this 
for a variety of t wo, three and four-quark operators. Our results are pre- 
sented in Table [^.IHHj There we present the renormahzation factors for each 
operator (continuum divided by lattice) as computed in perturbation the- 
ory, and in perturbation theory but with the (nonperturbative) mean-field 
factors removed. Perturbative expansions are much improved by ej^ii^cting 
the mean-field factors, particularly for operators with lots of fields 



(13) 



The results here are all for massless or nearly massless quarks. Tadpole- 
improved operators and actions for heavy quarks are discussed in [|l^, ^ , for 
Wilson quarks, and in [22|, for nonrelativistic quarks (NRQCD). 



5 Implications for Scaling 

A key issue in QCD concerns the onset of asymptotic or perturbative scaling: 
how small must the lattice spacing be before the variation of the coupling 
constant is perturbative. The variation of the coupling with changing lat- 
tice spacing is determined by the beta function, which, at short distances, is 
a perturbative quantity like any other. If perturbation theory successfully 
predicts a range of short-distance quantities, it is likely that it also correctly 
predicts the beta function. Thus our results in Section ^ provide indirect 

(12) 'pjjggg results are for Wilson parameter r = 1; very similar results arise for r = 1/2. 

(13) 'pjjg continuum operators used in this comparison were defined using the MS scheme. 
Our choice of normalization scale, fi = 1/a, was somewhat empirical; a more system- 
atic determination of the appropriate scale is possible using a variation of the techniques 
discussed in Section ^. Also, there is another obvious nonperturbative procedure for nor- 
malizing the operator for lattice quarks. The quark field's normalization should be roughly 
the square root of the normalization of either the vector or axial- vector current since these 
currents are conserved (or partially conserved). Inspection of our table indicates that us- 
ing the average normalization of the two currents to define the quark normalization gives 
even better results than those shown there. 
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operator 




perturbative 


tadpole-improved 




m 


(1 


— 1.39 ay) 


2kc 


(1 


- 0.03 ay) /4 




m 


(1 


- 2.40 ay) 


2kc 


(1 


- 1.03 ay) /4 




w 


(1 


-2.19ay) 


2kc 


(1 


- 0.82 ay) /4 




if 


(1 


— 1.68 av) 


2kc 


(1 


- 0.31 ay) /4 




ll 


(1 


- l.SOay) 


2f^c 


(1 


- 0.44 ay) /4 
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(1 


-2. 78 ay) 


(2kc)^-' 


(1 


- 0.73 ay) /8 






(1 


- 2.89 ay) 




(1 


- 0.84 ay) /8 






(1 


-2.76ay) 




(1 


- 0.71 ay) /8 




n 


(1 


— 3.81 av] 


{2Kcf 


(1 


- 1.08 ay) /16 






(1 


— 3.66 ay' 


{2Kcf 


(1 


- 0.93 ay) /16 




23] 


(1 


- 3.78 ay) 




(1 


- 1.05 ay) /16 


Tp'>p02lp1p 


H 


(1 


— 3.55 av] 


(2a€c)' 


(1 


- 0.83 ay) /16 



Table 3: Renormalization constants relating continuum (MS at = 1/a) and 
lattice operators for massless quarks. Results are for the ratio of continuum 
to lattice matrix elements, both without and with tadpole factors removed. 

evidence in support of perturbative scaling. Our results also test the scaling 
properties of the coupling constant directly. This is because at each /? we 
measure the coupling ay at q = 3.41/a, using data for the very ultraviolet- 
divergent plaquette, and then we perturbatively evolve ay down to scales 
ranging from 0.4/a to 2.8/a to compute estimates for a variety of less ultra- 
violet quantities. (Note that ay at /? = 6 increases by more than 50% when 
evolving from q = 3.41/a down to g = 1.09/a, the scale for Creutz ratio X22-) 
The success of our many perturbative estimates is compelling evidence that 
coupling-constant evolution is mostly perturbative for all /?'s down to 5.7, 
and possibly even for lower ones. Of course, this discussion only applies to 
the ay and a^^ definitions of the coupling; aiat is poorly behaved, but also 
largely irrelevant given our new perturbative techniques. 

The q-dependence of av{q) is readily extracted from our data. The 
results for three values of f3 are shown in Fig. |8|. To obtain these plots, 
we fit second-order expansions in av{q*) to Monte Carlo data for the six 
smallest Creutz ratios, and for the logarithms of the six smallest Wilson 
loops. The value of ay (g*) obtained from the fit for each quantity is plotted 
versus the q* for that quantity. The g*'s for the twelve quantites used here 
range from 0.43/a (for XAi) to 3.41/a (for — In Wn) — about a factor of eight. 
For comparison, we have included the (two-loop) perturbative prediction 
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for aviq) (solid line), arbitrarily normalizing ay so that the curve passes 
through the data point for —In 1^22-, The data are quite consistent with 



perturbative scaling, even at /3 = 5.7.1(1111 Note that statistical errors in the 



Monte Carlo data are negligible here; the fluctuations visible in the plots are 
due to uncalculated third-order terms in perturbation theory, which differ 
from quantity to quantity, and, at the lowest P, to nonperturbative effects. 
(The onset of the long-distance area-law in the logarithms of the Wilson 
loops is apparent in the plot for P = 5.7, although the effect is not all that 
large even for the 3x3 loop.) 

Our conclusion, that scaling is asymptotic even at /? = 5.7, contradicts 
standard lore. This lore derives from studies of scale invariant ratios of Aiat 
with physical quantities like the deconfining temperature Tc or the string 
tension a. Such ratios, which should become independent of /? at the onset 
of asymptotic scaling, show considerable variation with P for /3's less than 
6.2. This is illustrated by the upper plots in Fig. ^ These show the IP-IS' 
mass splitting AM divided by Aiat for the Tp and T [25| meson families, 



as well as -^/a/Aiat ||2^, for a range of /3's. Scaling violations of order 30-40% 
are readily apparent between /? = 5.7 to /? = 6.1. However, the situation 
changes completely if we replace the Aiat's in these ratios by Ay's determined 
(perturbatively) from the ay's we used in Section ^. When compared with 
Ay, the data are consistent with asymptotic scaling to within a few percent. 

Exact scaling is not expected for any physical quantity. There are cer- 
tainly finite-lattice-spacing errors in each of the measurements we use here. 
These errors have been analyzed carefully for the T data they result in 
roughly 10% scaling violation over the range shown. Errors for the other 
two quantities are probably smaller since T's are the smallest mesons. Note 
that the IP-IS splitting in quarkonium mesons is one of few hadronic mea- 
surements that is suitable for studying scaling. This is because the splitting 
is almost completely insensitive to the heavy-quark's mass, and so depends 
only upon the coupling constant. (This is also why the '0 results shown are 
nearly indistinguishable from the T results.) 

The evidence suggests that physically interesting quantities like mass 
splittings or the string tension scale perturbatively even at low /3's. The 



Similar results are reported for the SU2 lattice theory in although the finite- 
scaling technique used there is quite different from our procedure. That study probes 
different scales by examining a single quantity on a series of lattices with different lattice 
spacings. Our study probes different scales on single lattices by examining a variety of 
quantities, some more ultraviolet than others. In both cases the evolution of the coupling 
constant is tracked over a large range of scales. 
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Figure 8: Values of ay (g) for a range of g's as determined from lattice QCD 
measurements at various /?'s. The data points (circles) are measured values 
(with negligible statistical errors) obtained by fitting second-order perturba- 
tion theory to Monte Carlo simulation data for various short-distance quan- 
tities. The solid line shows the variation in ay(g) expected from two- loop 
perturbation theory. 
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Figure 9: Ratios of the square root of the string tension a and the IP-IS 
mass sphttings AM for -^'s and T's with Ai^t (top row) and with Ay (bottom 
row). 
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Figure 10: Ratios of the square root of the string tension a with the IP-lS* 
mass sphttings AM for ■^'s and T's. 



problem with the standard lore is that aiat does not scale perturbatively (at 
least through two- loop order). This is clear from our studies of perturbation 
theory. These studies also indicate that renormalized couplings like ay are 
perturbative, and this is why the ratios with Ay scale so well. Of course, 
ratios of physical quantities should scale properly as well, and they do (see 
Fig. 

It has been apparent for some time that the deconfining temperature 
scales better when analyzed in terms of a modified coupling constant sim- 



ilar to ours. 1 27 1 Now it is apparent that the modified coupling constant is 
just a continuum coupling constant like ay- Furthermore it is clear that the 
failure of scaling was intimately related to the lack of convergence of per- 
turbation theory for short-distance quantities like Kc or X22- Both problems 
are resolved by replacing a\^t with ay- 



6 Summary 

The use of lattice perturbation theory in conjunction with simulations has 
been hampered by two problems: 
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• expansions in powers of the bare lattice coupling aiat consistently un- 
derestimate perturbative effects, sometimes by factors of 2 or 4; 

• expansions for many quantities (and particularly renormalization con- 
stants for lattice operators) have large coefficients due to tadpole dia- 
grams and consequently converge poorly, if at all. 

We have addressed both problems in this paper. We have shown here that 
lattice perturbation theory works well when a proper coupling constant is 
used; and it can be made about as convergent as the continuum theory by 
systematically removing tadpole contributions. 

The first problem is remedied by replacing aiat with a renormalized cou- 
pling constant like av{q*), where scale q* is customized (in a predetermined 
way) to the quantity under study. The coupling constant ay is defined in 
terms a physical quantity, the heavy-quark potential, and it can either be 
measured (Section |^ or it can be determined from the bare lattice cou- 
pling aiat using formulas from mean-field theory (Section |4.2| ). Perturbation 
theory, when expressed in terms of ay{q*), is remarkably effective even at 
P = 5.7. 

The second problem, large tadpole-induced renormalizations, is reme- 
died by simple redefinitions of the basic operators used to define the lattice 
theory. Every in a naive lattice operator is replaced with Ufj,/uo, where uq 
is a measured constant representing the mean value of the link (Section Q); 
and every renormalized low-mass (Wilson) quark field t/2k^'0 is replaced 
by ip/2. The new operators obtained this way are rescaled versions of the 
naive operators. Their normalizations are very close to those of their con- 
tinuum analogues; renormalization constants for composite operators built 
from these tadpole-improved operators have perturbative expansions that 
are far more convergent. Tadpole improvement is essential for operators, 
like the cloverleaf operators for F^^, that involve many links; without it nor- 
malizations are wrong by as much as a factor of 2, and perturbation theory 
becomes useless. 

Our examples suggest that lowest-order perturbation theory in ay gives 
results for short-distance quantities that are typically correct to within 10- 
20% at /3 = 6. Expansions in aiat can be off by factors 2 or 4 at the same 
f3. Adding in higher-order corrections usually reduces errors by factors of 
2-5 for ay expansions, and by very little for aiat expansions. In many 
situations, ay expansions can be made still more accurate through tadpole- 
improvement, where powers of the mean-field parameter uq are factored out 
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of the expansion leaving behind a more convergent series. Our tadpole- 
improved one-loop formula (Eq. ^3|) for the critical value of the hopping 
parameter in Wilson's quark action, for example, is about as accurate as 
the best numerical determinations of this quantity. Finally, our procedure 
for determining the proper scale q* for the coupling consistently leads to 
excellent expansions, although expansions in ay(7r/a) for quantities defined 
over one or two lattice spacings usually give errors that are within a factor 
of 2-3 of those obtained with ay{q*). 

The fact that perturbation theory seems to be working at /? = 5.7 implies 
that asymptotic or perturbative scaling should also work. We verified this 
here, for a number of physical quantities, by comparing their dependence 
on (3 with that of the scale parameter Ky for the renormalized coupling ay ; 
scale invariant ratios of these quantities with Ky (as opposed to Aiat) showed 
little variation all the way down to (3 = 5.7. These results suggest that 
the lattice spacings used in current simulations are small enough for reli- 
able studies of QCD. Indeed, if anything, they are unnecessarily small. It 
is probably much more cost effective to simulate QCD at /? = 5.7, while 
removing the 0{a,a^) errors that are important by correcting the action. 
Previous efforts at improving lattice actions have not been too successful; 
but these relied upon the use of expansions in aiat , and naive lattice opera- 
tors. The perturbative quantities we examine in this paper are very similar 
in character to the new coupling constants that appear in corrected actions 
and operators. Thus our success in computing these quantities (to within a 
few percent in most cases) suggests that the use of ay pertrubation theory 
and tadpole-improved operators to correct the action will be much more 
successful. The potential savings in computer resources make it imperative 
that this possibility be thoroughly investigated. 

Appendix 

This appendix presents data for some of the figures in tabular form. 
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Perturbation Theory 


M. C. 




aiat "MS (9*) 


av{q*) 


Data 


5.7 


0.081 0.161 


0.191 


0.176 


6 


0.077 0.136 


0.157 


0.139 


6.4 


0.072 0.118 


0.133 


0.117 



Table 4: (1 — -gTiU) (Landau gauge) — the expectation value of the trace 
of a link in Landau gauge calculated in first-order perturbation theory and 
by Monte Carlo simulation for various /3's. Statistical errors in all data 
presented are of order one in the last digit quoted or smaller. 



p 


Perturbation 


Theory 


M. C. 




"lat 




av{q*) 


Data 


5.7 


-0.46 


-0.93 


-1.12 


-1.04 


6 


-0.43 


-0.78 


-0.91 


-0.80 


6.1 


-0.43 


-0.75 


-0.86 


-0.78 


6.3 


-0.41 


-0.70 


-0.79 


-0.70 



Table 5: anic — the critical quark mass rric for Wilson quarks with r=1.0, 
calculated in first-order perturbation theory and by Monte Carlo simulation 
for various /3's. Statistical errors in the simulation data are of order 2 in the 
last digit quoted. 



p 




First Order 




Second Order 


M. C. 




"lat 


"Ms 


ay 


"lat 


"ms 


ay 


Data 


5.7 


0.10160 


0.29790 


0.40100 


0.15429 


0.33181 


0.35347 


0.37343 


6 


0.09652 


0.23170 


0.28560 


0.14407 


0.25225 


0.26150 


0.26558 


6.2 


0.09341 


0.20840 


0.24990 


0.13794 


0.22497 


0.23146 


0.23317 


9 


0.06435 


0.09749 


0.10490 


0.08548 


0.10112 


0.10167 


0.10173 


12 


0.04826 


0.06399 


0.06703 


0.06015 


0.06556 


0.06570 


0.06574 


18 


0.03217 


0.03826 


0.03930 


0.03746 


0.03882 


0.03885 


0.03885 



Table 6: X22 — the expectation value of the Creutz ratio X22 calculated in 
first-order and second-order perturbation theory and by Monte Carlo sim- 
ulation for various /3's. Statistical errors in the Monte Carlo results are 
negligible. 
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n 


Perturbation Theory 


M. C. 




"lat 




av{q*) 


Data 


2 


0.13794 


0.22497 


0.23146 


0.23317 


3 


0.05207 


0.09467 


0.10312 


0.11348 


4 


0.02525 


0.06283 


0.07121 


0.06793 


5 


0.01512 


0.03881 


0.04780 


0.04949 



Table 7: Xnn — diagonal Creutz ratios Xnn as computed in second-order per- 
turbation theory and by Monte Carlo simulation at /3 = 6.2. Statistical 
errors in the Monte Carlo results are negligible. 



a q 


0.5 


0.8 aq* = im 


1.5 


2.0 


3.0 


X22 


0.174 


0.227 0.231 


0.228 


0.222 


0.212 



Table 8: Perturbative predictions for Creutz ratio X22 using expansion pa- 
rameter av{q) with various g's at /3 = 6.2. Monte Carlo simulation gives 
X22 = 0.233. 



a q 


1 


2 a q* = 2.65 


3 


4 


6 


- In W22 


1.074 


1.446 1.485 


1.492 


1.494 


1.474 



Table 9: Perturbative predictions for — In Tl 22 using expansion parame- 
ter av{q) with various g's at /? = 6.2. Monte Carlo simulation gives 
W22 = 1.527. 



p 


measured 


mean field 


perturbative 


5.7 


0.188 


0.168 


0.148 


6 


0.156 


0.145 


0.134 


6.2 


0.143 


0.135 


0.127 


9 


0.074 


0.073 


0.072 


12 


0.050 


0.050 


0.049 


18 


0.031 


0.030 


0.030 



Table 10: ayiTr/a) as measured (from — InT^n), as computed from the bare 
lattice coupling using nonperturbative mean-field theory, and as computed 
in perturbation theory. Statistical errors in the Monte Carlo results are 
negligible. 
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/5 


AM^/Alat 


AMx/Alat 


\/o^/Alat 


AM^/Av 


AMt /Ay 




5.7 


130(9) 


136(15) 


134(1) 


1.36(9) 


1.42(15) 


1.40(1) 


5.8 






120(1) 






1.40(1) 


5.9 


105(6) 




110(1) 


1.31(7) 




1.37(1) 


6.0 




105(4) 


105(1) 




1.36(5) 


1.38(1) 


6.1 


96(6) 






1.31(7) 







Tabic 11: Scale invariant ratios for the IP-IS mass differences AM for ip^s 
and T's and the square root of the string tension a. 





q* 


Fitted aviq*) 

f3 = 5.7 (3 = 6.2 (3 = 9.0 


-InWu 


3.41 


0.183 


0.140 


0.073 


-lnH^i2 


3.07 


0.196 


0.146 


0.074 


-InT^is 


3.01 


0.200 


0.147 


0.075 


- In W22 


2.65 


0.222 


0.155 


0.076 


- In W23 


2.56 


0.234 


0.159 


0.077 


- In W33 


2.46 


0.251 


0.163 


0.077 


X2.2 


1.09 


0.353 


0.208 


0.087 


X2.3 


1.10 


0.390 


0.215 


0.087 


X2.4 


0.97 


0.489 


0.235 


0.089 


X3.3 


1.07 


0.449 


0.227 


0.088 


X3.4 


0.80 


0.594 


0.262 


0.092 


X4.4 


0.43 


1.642 


0.370 


0.101 



Table 12: av{q*) as determined by fitting second-order perturbation theory 
to Monte Carlo simulation results for logarithms and Creutz ratios of small 
Wilson loops. Statistical errors in the Monte Carlo results are negligible. 
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